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Motivation 


Flight Feasible Trajectories will 
Model Realistic In-Flight Thermal States: 

• Allow for increased accuracy in Thermal Protection System sizing 
(potential mass savings) 

• Reduce the number of design cycles required to close an entry 
spacecraft design (potential cost savings) 


Novel Research Objective 

Develop a planetary guidance 
algorithm that is adaptable to: 

-Mission Profiles 

-Vehicle Shapes 

for integration into vehicle 
optimization. 


Altitude 
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De-Orbit 


Sample Concept of Spaceflight Operations 


* Adapted graphic from NASA Johnson Space Center 

De-Orbit 






Exploration: 

Vehicle completes mission 
over several day or weeks 


Launch to: 

• Earth Orbit 
Planetary Body 
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Planetary Entry Spacecraft Design (cont’ d) 


Mid - Low L/D Spacecraft 



* Orion Capsule 
www.nasa.gov 


a - variable bank angle 
a - fixed angle of attack 


ift % 



grav«y 


* MSL Capsule 
Prakash et al., 
NASA JPL 


High L/D Spacecraft 


Ellipsled 
Garcia et al., 
AIAA Conf. Paper 




a - variable bank angle 
a - variable angle of attack 






* NASP 

AIAA 2006-8013 


* HL-20 
AIAA 2006-239 
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Multi-Disciplinary Design, Analysis, and Optimization 

(MDAO) 


Computer Generated 
Spacecraft Models 


Available 

Descent 

Technologies 

Un/manned 

Planetary 

Models 

Mission 

Profile 




> 
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Minimize: 

Heat Rate (Trajectory/Shape) 
Ballistic Coefficient (Shape) 




Aerodynamic (C D , C L ) & 
Aerothermodynamic (q) 
Databases 


'Decoupled 
— ' Iterations 


Decoupled 

Iterations 


Flight Feasible 
Trajectory Database 

(replace Traj. Opt.) 




Guidance, 
Navigation, & 
Control 
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Trajectory Optimization vs. Guidance 



Trajectory Optimization 

Guidance 

Constraints 

Multiple included 

Minimal included 

Objective 

Any variable of interest 

Target specific 

Solution 

Purely numerical 

Combination of numerical and 
analytical 

Time to Solution 

Minutes to hours 

Seconds 

Guaranteed Solution 

No 

Must enforce that a solution is 
found 

Parameter Changes 

Handles large parameter 
changes 

Handles parameter changes 
that are relatively small 

Result 

Nominal Trajectory - not 
always realistic control 

Flight Feasible Trajectory 
with realistic controls 
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Guidance Development Trade-Offs 


Adaptability 

Numerical formulation for adaptability to different vehicles and missions 

without significant changes 


Rapid Trajectory Generation 
Analytical driving function keep time to a solution low 


Minimize Range Error & Heatload 
Optimal Control theory to introduce heat load as an additional objective 
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Guidance Development Criteria 

Guidance Specific (In-Flight) 

• Determine flight feasible control vectors (control rate/acceleration 
constraints) 

• Be highly robust to dispersions and perturbations 

• Include a minimal number of mission dependent guidance 
parameters 

Vehicle Design Specific 

• Be applicable to multiple mission scenarios and vehicle dispersions 

• Manage the entry heat load in addition to achieving a precision 
landing 


Types of Guidance Techniques 


Reference Tracking Only - follow a pre-defined track 

In-flight Reference Generation & Tracking - Generate a real-time 
reference trajectory and follow that track 

In-flight Controls Search - One dimensional search, usually solving 
equations of motion numerically 

In-flight Optimal Control - Requires numerical methods to meet some 
cost function 
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Types of Guidance Formulations 



Analytical Guidance 

Numerical Guidance 

Advantages 

• Simple to Implement 

• Computation time minimal 

• Solution Guaranteed 

• Accurate trajectory solutions 

• No simplifying assumptions 
(possibility of multiple entry cases 
to be simulated with few 
modifications) 

Disadvantages 

• Simplifications reduce accuracy 
of the trajectory solution 

• Formulation tied to a specific 
entry case 

• Convergence is not assured 

• Convergence is not timely 
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Novel Approach to Guidance for MDAO 

Real-Time Trajectory Generation and Tracking 

Adaptability 

Numerically solve entry equations of motion 
Use generalized analytical functions to represent the reference 



Rapid Trajectory Generation 

Use analytical driving function keep time to a solution low 
Use Single Opti mal Control Point with Blending 


Adaptation of Shuttle Entry 
Guidance Techniques 


Minimize ±w™ge Error & Heatload 
Optimal Control theory used to introduce heat load objective 



Adaptation of Energy State 
Approximation Techniques 
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Skip Entry Critical Points 

Test Case: Orion Capsule, L/D 0.4 Begin with 1 st Entry portion of the 
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Trajectory Simulation Validation 


{ Simulation of Rocket Trajectories (SORT) 
Developed by NASA Johnson Space Center for 
Space Shuttle Launch/Entry Simulations 
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Flight Dynamics 



b - body fixed coordinate 







ECF 


0 - longitude 
(|) - latitude 
y - flight path angle 
vj / - azimuth 


Earth Centered Fixed 
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Trajectory Modeling 


r = V sin y 

V cos y sin y 


V = -D - gsiny + Q 2 r cos ^(siny cos <j>- cosy sin ^ cos y/) 


0 = 


</>= 


r cos^ 

V cos y cos y/ 


Y = — 
V 


¥ = — 
V 


Lcosa + cosy 


f y2 

g 

r 


+ 2Q V cos (f> sin y/ + Q 2 r cos ^(cos y cos (j) + sin y cos y/ sin (f) 


Lsincr V 2 . . . . ,\ rQ 2 . . . . 

1 cos y sin y/ tan (p - 2Ll V (tan y cos y/ cos (p — sin (p) H sin y/ sin (p cos (p 

cos y r cos y 


State Variables 
r - radial distance 
V - relative velocity 
0 - longitude 
<|> - latitude 
y - flight path angle 
\\f - azimuth 


Control Variables 
a - bank angle 
a - angle of attack 


Vehicle and Planet Vanabhs 
L, D - Lift, Drag Acceleration 
g - gravity 

Q - Earth ‘s Rotation 
p - atmospheric density 
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General Entry Guidance Block Diagram 
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Control Solution: Shuttle Entry Guidance Adaptation 

Shuttle Entry Guidance (SEG) Concept: Temperature Phase 

• Reference Tracking Algorithm, Closed Form Solution 
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Control Solution: Shuttle Entry Guidance Adaptation 

Improvements on Shuttle Entry Guidance “Drag Based Approach” 

•Increase # of segments 
•Increase order of polynomial 
•Change Atmospheric Model representation 
•Modify flight path angle representation 

Challenges with Drag Based Approach 

• Discontinuities between segments 

• Increasing # of coefficients for storage with increasing segments and/or 
order 

• Effect of small flight path angle assumption unknown 

• Formulations are derived from 2DOF Fongitudinal EOMs 
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Control Module: Shuttle Entry Guidance Adaptation 

Sensitivity to atmospheric non-linearity is significant during initial and final 
segments. Need an Alternative Analytical Eauation! 



20 


Automated Selection of Transition Events 

Framework: 

- Allows for adaptability 

- Automated generation of Reference Trajectory 

- Open loop 


Study Objective: Define bank profile for trajectory phases 


Phase 

Bank Description 

Entry Interface to 
Guidance Start 

Constant Bank 

Guidance Start to 
Guidance End 

Trajectory Solver 

Guidance End to 

Exit 

Linear Transition to Meet 2 nd Entry Bank 

Exit to 2 nd Entry 

Attitude Hold 
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Automated Selection of Transition Events 


• Metric to determine best trajectory: lowest range error, lowest heat load from 
El to 2 nd Entry, and bank transitions 
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Automated Selection of Transition Events 

Study Results: 


Phase 

Bank Description 

Entry Interface to Guidance Start 

Constant Bank = 57.95° 

Guidance Start to Guidance End 

Trajectory Solver 

{0.12 0.11} G’s 

Guidance End to Exit 

Linear Transition to Meet 2 nd Entry Bank 
Linear Transition Velocity: 23,784.65 ft/s 

Exit to 2 nd Entry 

Bank Attitude Hold = 70° 
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General Entry Guidance Block Diagram 
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Targeting Algorithm Development 

When is Targeting Activated? 


1 .Overshoot - Vehicle is predicted to fly wav past target 
2. Undershoot - Vehicle is predicted to fly short of the target 

How to find a set of controls to Correct Over/Underhoot? 

Adapt Energy State Approximation Methods’. 

Optimal control method that replaces altitude U r 2 

and velocity with specific energy height (e) 2 g 

Advantages: Allows for a compact set of analytical equations 

Add heat load to the range error objective function 

Disadvantage: Optimal control formulations may not converge to a solution 

Solution: Derive a localized optimal control point instead and blend 
back reference trajectory 
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Targeting Algorithm Development 

Must Relate Euler-Lagrange Equation 

A = = tan a* cos 7 

A7 < 0 

To Reference Trajectory Variables 


L_ 

D 


COS cr = - 

total 


ref 



( v 2 ^ 


Vj ref - cos y 

y r 

g 

~C 7 (y) 


l ^ J 



Using trigonometry and other manipulations, the control equation is 
found 


L_ 

D 


total 


1+ (^) 


^7rej~ cos 7^~~~5 f ^ —201/ cos sin P 2 r cos 0(cos7cos tft+sin 7 cos iff sin 


D 


aprx 
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Targeting Algorithm Development 



Realtive Velocity [ 
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Targeting Algorithm Development 

Targeting Technique 1 - Design Space Interrogation 
C* - drag/density ratio coefficient 
d\ - change in Lagrange multiplier 
dV - change in relative velocity at next point 

Targeting Technique 2 - Design Space Interrogation 
d\ - change in Lagrange multiplier 

dV i - change in relative velocity halfway to curve fit end point 
A (dE) - second order change in energy 
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D/p |ft 3 /slg • fVs 2 )| 


Targeting Algorithm Development 


Targeting Technique 1 - Design Space Interrogation 



dk |N1)J 


Targeting Algorithm Development 


FPA Dispersion - Undershoot 



Undershoot Case - Open Loop % Range Error 


dV [lVs| 



Undershoot Case - Closed Loop Range Error |km] 


dV | ft/s 1 


Undershoot Case - I leal load Contour |J/cm 2 | 



Icinl] xp 


Targeting Algorithm Development 





s 

c ^ 
XI 


r> i 

T3 


dV iU/s| 


FPA Dispersion - Overshoot 

Overshoot Case - Open Loop % Range Rnor 


Overshoot Case - Closed Loop Range Error [km] 


dV (ft/s | 


Overshoot Case - I lea I load Contour [J/cm 2 ] 
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Icin] xp 


Targeting Algorithm Development 

Aerodynamic Dispersion - Overshoot 

Overshoot Case - Open Loop Range Error [%| Overshoot Case - Closed Loop Range Error |km| 




Overshoot Case - Ileatload Contour |J/cm 2 ] 






Shape Optimization Analog 


Computer Generated 
Spacecraft Models 



Geometry #1: C L = 1.70, C D = 3.4 
* Geometry #2: C L = 1.90, C D = 3.8 



Geometry #3: C L = 1.95, C D = 3.9 

Current Guidance Algorithms - Robust to 
~20% aerodynamic dispersions 



Must exceed 20% to demonstrate 
potential for integration into MDAO 

ANALOG: Changing 
angle of attack 
disperses C L and C D 




Range Miss [nmi] 


Targeting Algorithm Development 

Guidance Algorithm for Comparison - Apollo Derived Final Phase Guidance 
Reference Tracking to a stored trajectory database, function of relative velocity 


Performance Results - Threshold Miss Distance, 1 nmi 


800 


600- 


400- 


200 - 


43.62- 


7.819k 


a = 167° 


1 


A a=l° 


l 


a =158° = 0 

Range Error a - 153 

Starts to Grow Reference Case 


. . I 


•A %, 



0.15 


0.2 


0.25 0.3 

L/D total P®] 


0.35 


0.4 


0.45 


34 


Targeting Algorithm Development 


Targeting Technique 1 - Targeting Procedure 

1 . Guess a value for dk 

2. Iterate on dV using secant method to converge on a zero 
range error trajectory 

3. If no solution is found, dk is incremented and the iteration 
repeated 

4. Solution is then flown in flight simulation 


Bank Angle [deg] 


Targeting Algorithm Development 

Targeting Implementation, 1 st and 2 nd Phase - Results 

a =152°, L/D -0.418 a =153°, L/D -0.402 a -154°, L/D -0.386 



time [s] 


Targeting Algorithm Development 

Targeting Technique 2 

V r 

Use Energy Height e - - 1 — + h to determine Control Point \V new , $ new ] 

2 go 

Undershoot — > energy dissipating (de/dt) too fast 
Overshoot — ► energy dissipating (de/dt) too slow 

Since Velocity is an independent variable 
and a pseudo control de/dV is examined 




Targeting Algorithm Development 

Targeting Technique 2 

Recall the equation for the ratio of drag acceleration to density: 


D 

P 


-Extract altitude and velocity from [dVi, A (dE)] to find O 


new 


st 

, 8 Notional - Undershoot Case for 1 Entry Phase 
x 10 



Targeting Algorithm Development 

Targeting Technique 2 - Design Space Interrogation 


d A 

Lower Limit 
0 

Upper Limit 

* 

Incr. 

■ ■■■■■■ 

units 

uuuuuuuu 

ND 

dV l 

0 

1524 

Predict 

m/s 

A (dE) 

0 


Predict 

m 


^minfmax — tail <7 min /max COS 

Limit are trajectory dependent 

and control system dependent _ /* _t \ 

delimit — “F y^min/max ^ old ) 


Dispersion Cases: ^ ^ 

L/D Dispersion 

Target Miss 

1 st Phase Only Nominal 

0.4 (0%) 


152 

0.42 (+5%) 

Undershoot 

162 

0.28 (-30%) 

Overshoot 

165 

0.23 (-43%) 

Overshoot 

167 

0.2 (-50%) 

Undershoot 
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Targeting Algorithm Development 

Design Space Interrogation, Results: Range Error [%] 


1A 4 a = 1 52°, Undershoot 

x io 




Targeting Algorithm Development 

Design Space Interrogation, Results: Heatload [J/cm A 2] 


x io 4 a = 152°, Undershoot 



dV [m/s] 



a = 165°, Overshoot 

X 10 



a = 167°, Undershoot 

x 10 



dV [m/s] 


Targeting Algorithm Development 

Design Space Interrogation, Results: Bank Rate [deg/s] 


a = 152°, Undershoot 

x 10 




10 a = 162°, Overshoot 



300 350 

dV | m/s | 


x 10 -» a = 167°, Undershoot 
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Targeting Algorithm Development Results 


Dispersions - 

Apollo Derived Guidance = -20% dispersion 
MDAO Algorithm = -43% dispersion 

Managing heatload may be a challenge for dispersions greater than 20% 
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Conclaamasi Specific (in-Fiight) 

V Determine flight feasible control 



vectors (control rate/acceleration 
constraints) 

o 

Be highly robust to dispersions 
and perturbations 

V 

Include a minimal number of 
jraissj^de^ndent guidance 
/p^jameters ref 

Vehicle Design Specific 

• 

Be applicable to multiple 

O/V mission scenarios 
V vehicle dispersions 


o 

Manage the entry heat load in 
addition to achieving a precision 
landing 


Acknowledgemets 

University of California, Davis 

Dr. Nesrin Sarigul-Klijn 
Dissertation Chair 

Dr. Dean Karnopp 

Dissertation Committee Member 

UC Davis Mechanical and 
Aerospace Engineering Faculty 

UC Davis Mechanical and 
Aerospace Engineering Staff 


NASA Ames Research Center 

Dr. Dave Kinney 
Dissertation Committee Member 

Mary Livingston 
Supervisor 

Colleagues in Systems Analysis 
Office 


Thank You !!! 


45 


Questions? 
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Overview 


Background & Motivation 


MAPGUID Development 


Aerothermal Management 
During Guidance 


Guidance/COBRA 

Integration 


Closing Remarks 


i Elements of Spacecraft Design 

■> 

| Introduction to Planetary Entry Guidance 

j Dissertation Research Plan and Status 

A 

j MAPGUID Proposed Approach 
! Key Results #1 



■ Key Results #2 

■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■I 

i Key Results #3 



j Key Results #4 



j Proposed Approach 

j Key Results #1 



i Key Results #2 

j Proposed Approach 


■ Key Results #1 

■p 

■ Key Results #2 
j Key Results #3 

i 

■ Dissertation Findings and Status 


Big Picture: 

Spacecraft Design Process 
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MDAO Literature Review 


Vehicle Optimization and TPS Sizing 

Example Objective Function: qconv = 1.83 * 10 4 v@^(l - hw/Hsjffi 


Results 
• Most studi 
maximum 


What is Flight Feasible? 

• Reaches Target @ Landing Speeds 

• Control does not exceed system limits 


corresponding to 


• Used _ — — — - t rate 

Some studies use new trajectories, but there is no accounting for bank constraints or 
target accuracy 

None of these studies incorporated flight feasible trajectories 
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Proposed Approach to MDAO for Spacecraft Design 


Computer Generated 
Spacecraft Models 


Available 
Descent — 
Technologies 

Un/manned — 

Planetary 

Models 

Mission 

Profile 


■> 

-> 

■> 


Vehicle A Planetary 
Optimization [ ) Entry 

Guidance 


[up! 


Thermal 
Protection 
System 
(TPS) Sizing 


Structures 


1 


Aerodynamic (C d/ C L ) & 
Aerothermodynamic (Cf) 
Databases 

1 


I 

j Reduced 
Decoupled 
* Iterations 


Reduced 
I Decoupled 
j Iterations 


Flight Feasible 
Trajectory Database 


Guidance, 
Navigation, & 
Control 


51 


Trajectory Modeling for Design vs. 
In-Flight Trajectory Modeling 
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Planetary Entry Guidance Literature Review 


• High L/D, Earth: Space Shuttle, X-33, X40A 

• Most Robust: In Flight Trajectory Shaping with Reference 
Tracking 

• Least Robust: Reference Tracking Only 

• Low L/D, Earth: Apollo, Orion 

• Most Robust: In-Flight Controls Search 

• Least Robust: Reference Tracking Only 

• Other Planetary Entry Vehicles: MSR, MSL, Biconic 

• Flight Tested algorithms preferred 
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Planetary Entry Guidance Literature Review (cont’ d) 

Key Results 


Modem guidance algorithms: optimal control is potential framework, but 


convergence still an issue 
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Trajectory Optimization Literature Review 


Trajectory Optimization 

Traj - Nonlinear constrained optimization 
Mission - Sequential Quadratic Programming 
Energy State Method - Reduced Order Modeling, one dimensional 
parameter search 

Pseudospectral Methods - Combination indirect and direct 
method, mapping and discretization of domain 
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Trajectory Optimization Literature Review (cont’ d) 

Key Results 


Siridai^efgmid^^ itiaale increases with dimensionality 
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Introduction to Planetary 
Entry Guidance 
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Guidance Development Process 


Trajectory Design 


Trajectory Optimization 

-classical optimization 
-genetic optimization 

Bank Profile Scan 

Determine constraint boundaries to 
identify flight corridors 


Constant Bank 

Modulated Bank 

linear ramp, user-defined 
scans at each time step 


Reference Trajectory Found 


Guidance 

Development 


Identify Target: 


Trajectory Modeling in Real-Time 

-Orbit 

► 

to Predict Target Acquisition: 

-Landing Site 


-Profiles Stored in Tables 

-Descent Device 


-Analytical Equations 

Activation State 


-Numerical Integration of EOMs 

-Multiple 


-Hybrid Approach 


v_ 


Targeting Algorithm: 

-Parameterization 
-One-Dimensional Bank Search 
-Optimal Control Formulation 
-Analytical Equation Blending 
-Search Family of Trajectories 
-Hybrid Approach 


*Guidance must be robust to many dispersions: 
(Atmospheric properties, Aerodynamics properties, 
Navigational Inputs, Entry Interface Conditions, Mass, 
Control System performance, and many others) 


T 
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Baseline Vehicle & Mission 
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Case Study Parameters 


Vehicle 

Orion Capsule, L/D = 0.4 

Trajectory 

Skip Entry for Lunar Return 

Control 

Bank Angle only 

Atmospheric Model 

1976 Standard Atmosphere 

Gravity Model 

Central Force + Zonal Harmonics 

Aerodynamics 

C L , C D corresponding to Mach # 

Trajectory Simulation 

CBAERO Databases, function of 
Mach #, Dynamic Pressure, and 
Angle of Attack 

MATLAB Simulation validated 
against SORT Trajectories 


Trajectory Simulations Developed 

Open Loop Numerical Predictor- Corrector (NPC) Simulation 

Used to test guidance formulations 


3DOF Rotating Spherical Planet 


r = V sin y 


6 = 


V cosy sin y/ 
r cos(f) 




V cosy cos y/ 


V = —D - g sin y + Q 2 r cos ^(sin y cos (j> - cos y sin (j> cos y/) 


7 = 
Y = 


J_ 

V 

J_ 

V 


L cos cr + cosy 
Z sin cr V 2 


( y2 


\ 


g 


V 


+ 2Q.V cos^sin^ + Q rcos^(cosycos^ + sinycos^sin 


rQ 2 


cosy 


cos y sin y/ tan (f)-2Q.V (tan y cos y/ cos <f> - sin (f ) h sin y/ sin (f) cos <f> 


cosy 


Flight Simulation - Closed Loop Guidance Testing 

Using equations derived from Newton’ s 2 nd Law, dynamics of relative 
motion, and Earth Centered Inertial (ECI) coordinate system 
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Trajectory Solver Development 
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Control Solution: Shuttle Entry Guidance Adaptation 


Drag Curve Fit Accuracy 

# of stored 


Segments Order 


coefficients 


7(3) 

7(5) 

14 

7 


Irrational 

Irrational 

5 

2 


168 

105 

84 

21 



= C 2 V X2 +CjV Xl +C 0 V X ° 


¥ 




o Nominal 

* 10 pts/seg; 5th Order 
+ 10 pts/seg; 2nd Order 

* 20 pts/seg; 2nd Order 



X 

| 9 

X 


1 

1 1 

i i i i 

X 


6 1 „ , 1 1 1 — i -i— : 1 

3.51 3.5105 3.51! 3.5115 3.512 3.5125 3.513 

V(^) xl0 4 
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Control Solution: Shuttle Entry Guidance Adaptation 


Would Cubic Spline Interpolation work? h 


1 dP 1 dT\ _1 
P~dh~Tdh) 
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Targeting Algorithm Development 
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Altitude [ft] 


Targeting Algorithm Development 

Targeting Technique 1 - Trajectory Behavior to Full Set of Aerodynamic Dispersion 


Can Technique 1 find a trajectory that points toward correcting 
the range error? 



General Conclusions 
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Targeting Algorithm Development 


Pontryagin’ s Principle in Optimal Control 

Find Optimal Control u* an(l v *0 

for dynamic system x = fix, u , t) 

The optimal control V . /^ eo . raints including the Euler- 
Lagrange Equation: e = 


a// 


tl=U' + 



(!) 


= Fcos 7 A ' ■ u 


tan a* cos 7 


7 = 


X 

1 / 


A7 < 0 

L cos a -t- cos 7 -r zi ' v cos (p sin 0 

-|-f) 2 r cos 0(cos 7 cos 0 + sin 7 cos 0 sin 0)] 


$ 



L sin (7 
cos 7 


+ 



cos 7 sin 0 tan 0— 


2 QV (tan 7 cos 0 cos 0 — sin 0)+ 


rO 2 
cos 7 


sin 0 sin 0 cos 0 


Targeting Algorithm Development 

Targeting Technique 1 

\ t Determines new bank 

Kew = Add angT6t«\«urrent time step 

7 = / L cos cr + cos 7 ( V- — g J + 2 QV cos (p sin i[) 
— QV cos 0 ( cos 7 cos 0 + sin 7 cos 0 sin ^ 




rieiL’j&o-und — 


r-=l 

Jl+( cq ^ ) Vi ref —cos 7 ^ — 2 £lV cos <f> sin ip—i 1 2 r cos $(cos 7 cos 0+sin 7 cos ^ sin $) 


D — 

y D total 




old 


Calibrated for Each 
Dis^ersbd'Case 


t r _ 1 / , n ( y{^7r\ Determines Blended 

* initial — * current “r U.UM47 . , , ,, 

Trajectory that nulls 

^ re/,/ — h current (1 0.01 )tiHaT"IQe error 
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Targeting Algorithm Development 

Targeting Technique 1 - Design Space Interrogation 


• The blending technique exhibits potential to find new bank 
profiles that null the range error 

• The design space is constrained by control system limitations 

• There is a zero range error solution for each change in dX 
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Altitude [ft] 


Targeting Algorithm Development 

Targeting Technique 1 - Trajectory Behavior to Full Set of Aerodynamic Dispersion 



°o 200 40 The reference bank profile over-corrects with respect to the -&> 

dispersion of L/D 
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Targeting Algorithm Development 

Targeting Technique 2 

Now that the blended function is fully defined $bind — Bh 2 V 2 + Bb\V + Bb Q 


The following equation can be used to solve for: 


Hi, new — y 


L_ 

D 


total, i , I -rr^ P&Mntt + cos & sin 


cos 7 i ^ - g^j + QV* cos 0*(cos 7* cos + sin 7* cos ^ sin <p t ) 


The FPA rate table is shifted accordingly 
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Targeting Algorithm Development 

Design Space Interrogation, Results: Bank Acceleration [deg/s A 2] 


in 4 a = 152°, Undershoot 

x 10 7 



dV [m/s] 


io 1 a = 162°, Overshoot 



dV [m/s] 



x 10 ^ a = 167°, Undershoot 



dV [m/s] 
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Targeting Algorithm Development 

Targeting Technique 1 - Targeting Implementation, 1 st and 2 nd Phase 

1 . Guess a value for dk 

2. Iterate on dV using secant method to converge on a zero 
range error trajectory 

3. If no solution is found dk is incremented and the iteration is 
repeated 

4. Solution is then flown in flight simulation 


Performance Metric - 

Compare range of aerodynamic dispersions this algorithm can 
handle to the range of aerodynamic dispersions a heritage 
algorithm can handle. 
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Trajectory Solver Research Questions 


Can a simplification in the equations of motion be made without loss of 
accuracy? 

Can a simplification on flight path angle be made without loss of accuracy? 
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Simplified Equations of Motion Study 


3DOF Rotating , Spherical Earth (3RSP) 


¥ = 


J_ 

V 


Z sin cr V 2 . . rQ 2 

1 cos y sin y/ tan (p - 212 K (tan y cos y/ cos <p - sin (p) -\ 

cosy r j cos y 


3DOF Non-Rotating Spherical Planet 


¥ = 


V 


Z sin cr V 2 . . 

1 cos y sin y/ tan (p 

cos y r 


3DOF Non-Rotating Flat Planet 


¥ = 


V 


L sin cr 
cosy 


2DOF Longitudinal Equations (2LON) ^ 

h = Vsmy 
s = V cos y 
V = — D — g sin y 

. 1 


r = 


v 


Z cos <r + cosy 


( y2 




g 


V 


7 


V. Apollo and Shuttle 
Entry Guidance 


J 


sin y/ sin $ cos $ 


v 

Coriolis and 
Centripetal 

Acceleration 
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Altitude [ft] 


Simplified Equations of Motion Study (cont’ d) 
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Altitude [ft] 


Simplified Equations of Motion Study (cont’ d) 



Total Range [nmi] 
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Trajectory Solver Research Questions 


Can a simplification in the EOMs be made without loss of accuracy? 

Not for a skip trajectory 

Can a simplification on flight path angle be made without loss of accuracy? 
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Control Solution: Shuttle Entry Guidance Adaptation 


Kef = ~K 


A./ 2F c 

Z) . V C 


D 


ref 


D 


h s =RT/g Four Drag Segments (-30 points/segment) Quadratic Curve Fit for Drag 
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Altitude Rate [ft/s] 


Control Solution: Shuttle Entry Guidance Adaptation 
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Control Solution: Shuttle Entry Guidance Adaptation 

P = pRT 
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Control Solution: Shuttle Entry Guidance Adaptation 


Need to Resolve 1 st Segment to Capture Atmospheric Non-Linearity 

Tl 

IDEA: Curve fit drag with Mach Number jj ,=Z C t Ma l 


i=l 


1/H s = (dP/dh * P' 1 ) - (dT/dh * T' 1 ) 9 Segments (~1 5 pnts/seg) 6 th Order Drag Curve Fit 
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Control Solution: Shuttle Entry Guidance Adaptation 

Check Altitude Acceleration Approximation 

Segment 6 Segment 5 Segment 4 Segment 3 Segment 2 
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Trajectory Solver Research Questions 


Can a simplification in the EOMs be made without loss of accuracy? 

Not for a skip trajectory 

Can a simplification on flight path angle be made without loss of accuracy? 
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Range Prediction Sensitivity to Flight Path Angle 
Assumption 


• Apollo and Shuttle Entry 

guidance formulations 

approximate flight path 
angle (FPA) to be small: 


y « 1 rad 


and/ or y « 1 rad / s 


Why does this matter? 

• If predicted range does not equal the range to landing site then targeting is 
erroneously active 

• Are model reductions in the Trajectory Module and Control Module valid 
based on the nominal case? 
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Range Prediction Sensitivity to Flight Path Angle 
Assumption 


Case Studies: 

A. Apply y « 1 rad to Trajectory 
Module only 


B. Apply y « 1 rad to Controls 
Module only 


Trajectory Module 

NPC Solves 3DOF EOMs 

I I 


Controls Module 

I Drag and FPA Rate 

I Reference Trajectories 

L 


1 

I 

I 

I 

I 


C. Apply y « 1 rad Is to bank 
equation only 


L_ 

~D 


v,ref 


/XD 


ref _ 


Kfref~ COS r 


v 2 

r 


V 


\ 

g ~Cy( y) 

J 
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Range Prediction Sensitivity to Flight Path Angle 
Assumption 


Nominal 661.73 [nmi] 


Case 

Total Range [nmil j 

% Range Error 

T 1 

! Termination ! 

i a 

■ 

662.39 ! 

0.099% 

j Drag Limit j 

r 

B 

|_ 

649.74 j 

1.813% 

j Drag Limit j 

[c_ 

632.13 i 

4.474% 

j Velocity Limit \ 


Conclusion FPA approximation can be applied to the trajectory 

module, but not to the control module 
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